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ABSTRACT 

The Solar system beyond Neptune is believed to house a population of small 
primordial bodies left over from the planet formation process. The region up to 
heliocentric distance ~ 50 AU (a.k.a. the Kuiper Belt) may be the source of the 
observed short period comets. In this region, the phase space structure near orbital 
resonances with Neptune is of special interest for the long term stability of orbits. 
There is reason to believe that a significant fraction (perhaps most) of the Kuiper Belt 
objects reside preferentially in these resonance locations. This paper describes the 
dynamics of small objects near the major orbital resonances with Neptune. Estimates 
of the widths of stable resonance zones as well as the properties of resonant orbits are 
obtained from the circular, planar restricted three-body model. Although this model 
does not contain the full complexity of the long term orbital dynamics of Kuiper Belt 
objects subject to the full N-body perturbations of all the planets, it does provide 
a baseline for the phase space structure and properties of resonant orbits in the 
trans-Neptunian Solar system. 



1. Introduction 



It has long been conjectured, on both theoretical and observational grounds, that the 
outermost regions of the Solar system may contain a large population of small icy bodies. For 
example, on the basis of theoretical considerations of the genesis of the planetary system from the 
primordial Solar Nebula, Kuiper (1951) suggested that a remnant population of icy planetesimals 
left over from the planet-building process may exist at the current epoch beyond the orbit of 
Neptune. Whipple (1964)| and Bailey (1983) speculated on a massive comet belt as the source of 
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unexplained perturbations of Neptune's orbit [although this argument must now be discarded as 
the post- Voyager revisions in the planetary ephemeris no longer show any unexplained residuals in 
Neptune's motion (Standish (1993))]. Hamid et al. (1968) analyzed the orbital plane perturbations 
of comet P/Halley and concluded that any comet belt between 40 AU and 50 AU has a total mass 
less than 1M®. This limit still allows for very large numbers — perhaps 0(1O 9 ) — of cometary 
bodies. 

More recently, it has been suggested that the observed short period comets with orbital 
periods < 20 yr — the "Jupiter family" comets — originate in a belt of low-inclination bodies 
just beyond the orbit of Neptune, between 35 AU and 50 AU [ Fernandez (1980) , Fernandez fc Ip| 
(1983) 1- The older hypothesis that short period comets originate in a population of near-parabolic 



Oort Cloud comets which are perturbed into shorter orbits by the giant planets appears unlikely: 
Duncan et al. (1988)| , |Q uinn et al. (1990)] have shown that the orbital element distribution of the 



observed short period comets is inconsistent with a source in the nearly isotropic Oort Cloud but 
is compatible with a disk- like source in a trans-Neptune comet belt, which they call the "Kuiper 
Belt" . A possible member of the Kuiper Belt was first discovered in 1992 at a distance of 41 AU 
from the Sun [1992 QBi, reported in jjewitt &; Luu (1993) ], and several additional discoveries have 
been reported since suggesting a potential population of ~ 35,000 bodies larger than ~ 100 km 



[ Jewitt fc Luu (1995)[] . 



The dynamical structure of this population has been the subject of several recent theoretical 



studies [ Levison &: Duncan (1993) , Holman Wisdom (1993) , Malhotra (1995) |. The first two of 



these [ Levison &; Duncan (1993) and Holman fc Wisdom (1993) |, studied the long term stability 
of test particles in low-eccentricity and low-inclination orbits beyond Neptune, subject only to the 
gravitational perturbations of the giant planets in their present orbital configuration. They found 
orbital instability on short timescales (< 10 7 yr) interior to 33-34 AU, an intricate structure of 
interspersed regions of stability and instability in the semimajor axis range of 34 AU to 43 AU, 
and substantially stable orbits beyond 43 AU. The intricacy of the dynamical structure appears 
to be particularly acute near the locations of orbital resonances with Neptune. 

Any primordial trans-Neptune population of planetesimals was undoubtedly subject to 
dynamical evolution during the planet formation process, and the initial conditions assumed in 
the above studies are not necessarily representative of the state of the Kuiper Belt at the end of 



planet formation, as acknowledged in Holman fc Wisdom (1993) . The nature of the "dynamical 
sculpting" of the Kuiper Belt that would have occurred during the late stages of planet formation 



was the subject of a study by |Malhotra (1995)[ This study concludes that the giant planets' orbits 
would have evolved significantly and the Kuiper Belt would have been sculpted into a highly 
non-uniform distribution early in Solar system history: most of the primordial small bodies in the 
region beyond Neptune's orbit and up to approximately 50 AU heliocentric distance would have 
been swept into narrow regions of orbital resonances with Neptune, particularly the 3:2 and the 
2:1 orbital resonances which are located at semimajor axes of approximately 39.4 AU and 47.8 
AU, respectively. The orbital inclinations of most of these objects would remain low (typically 
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< 10°), but the eccentricities would be excited to significant values, typically 0.1 — 0.3. This 
structure would be largely preserved to the present epoch. 

Pluto (and its satellite Charon) have long been known to reside in the 3:2 Neptune resonance; 
the resonance libration protects this pair from close encounters with Neptune, despite their 
Neptune-crossing orbit (with eccentricity ~ 0.25). Several of the newly-discovered Kuiper Belt 
objects are also likely residents of this resonance | |Marsden (1995)| . (However, note that the orbital 
parameters of most of these objects remain rather poorly constrained pending future follow-up 
observations.) 

In any case, it is clear that the locations of Neptune's orbital resonances are particularly 
interesting with regard to the long term dynamics of the Kuiper Belt. 

In the region between Neptune's orbit and approximately 50 AU heliocentric distance, by 
far the most important orbital perturbations are due to Neptune alone, although the N-body 
effects of all the planets - particularly secular resonance effects - also surely play a significant 
role over the age of the Solar system. In this paper, I use the simplest possible dynamical model 
- namely, the circular, planar restricted three-body model - to determine the basic phase space 
structure near the locations of Neptune's orbital resonances exterior to Neptune's orbit. This 
model does not reflect the complexities that arise with (i) non-zero orbital inclination of the 
test particle, (ii) a realistic representation of Neptune's motion (its time- varying orbit), and (iii) 
the perturbations of the other planets. Some of these unmodeled effects can lead to important 
variations in the dynamics. In particular, as mentioned above, secular resonances in certain 
regions of phase space can drive the inclinations and eccentricities to large amplitudes and further 
complicate the dynamics near mean motion resonances. Nevertheless, this model is a reasonable 
first approximation because (a) the evidence from short period comets suggests a low-inclination 
source population, (b) in the trans-Neptunian space, there are no first order orbital resonances 
with any other giant planets, and (c) Neptune's eccentricity and inclination do not exceed ~ 0.025 
and ~ 2.5°, respectively, over billion year timescales under the effects of the mutual gravitational 
perturbations of the planets. The advantage of this simple model is that it is tractable: it 
allows for the calculation of 2-dimensional surfaces-of-section in which the phase space structure 

- in particular, the properties of resonance regions — can be easily visualized. This provides 
for relatively straightforward estimation of the widths of resonances and the determination of 
several interesting dynamical properties of resonant orbits, such as the libration periods and their 
dependence on orbital eccentricity and libration amplitude. To a considerable degree, most of 
these properties are preserved in the more realistic situation. In short, this model provides a 
reasonable "baseline" for the phase space structure in the trans-Neptune region. 

The rest of this paper is organized as follows. Section 2 describes the technical details of the 
planar, circular restricted 3-body model for the Sun, Neptune and a Kuiper Belt object, including 
a description of the information that may be gleaned from surfaces-of-section. Section 3 provides 
a detailed look at the phase space structure near the major Neptune resonances of interest for the 
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Kuiper Belt. Section 4 summarizes and discusses the results. 



2. The planar circular restricted three-body model 

The classical planar circular restricted three-body problem is a particular case of the general 
gravitational three-body problem of masses mi,m2,m3 defined by the following restrictions: 

1. the motion of all three bodies takes place in a common plane; 

2. the third body, 777.3, has zero mass; therefore it does not influence the motion of mi and 777,2; 
and 

3. the masses mi and 7772 describe circular orbits about their common center of mass. 

In the context of this paper, mi represents the Sun and 7772 represents Neptune. The system is 
made non-dimensional by the following choice of units: 

- the unit of mass is taken to be m\ + 7772; 

- the unit of length is chosen to be, ajy, the constant separation between 7771 and 7772 (i.e. the 
mean heliocentric distance of Neptune) ; 

- the unit of time is chosen such that the orbital period of mi and 7772 about their center of 
mass is 2ir. 

Then the universal constant of gravitation, G = 1, and the masses of the Sun and Neptune are 
1 — /i and /j, respectively, where /j = 7772/(7771 + 7772). For the Sun-Neptune system, we have 
fi = 5.146 x 1CT 5 , and this value is used throughout this paper. 

In a reference frame with axes (X, Y) rotating with mi and 7772 and with origin at their 
center-of-mass, the Sun and Neptune have fixed coordinates, (— /i, 0) and (1 — /i,0), respectively, 
and the third (massless) body has the following equations of motion: 

r l r 2 

(1) 

Y = -2X + Y-(1-^-^, 
where r\ and T2 are the test particle's distance to the Sun and to Neptune, respectively: 

n = [ ( x + /J ) 2 + Y 2 ] 1/2 , 

(2) 

r 2 = [(x-i + / x) 2 + y 2 ] 1/2 - 
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These equations of motion admit a constant of the motion, the Jacobi integral, given by 



C = X 2 + Y 2 -X 2 -Y 2 + 2(^ + ^1. (3) 

V n r 2 J 

No other constant of the motion is known. For future reference, we note that the Jacobi integral 
can be expressed in terms of the test particle's semimajor axis, a, and eccentricity, e: 



C = - + 2yXl-e 2 )+e>(/x). (4) 



Surfaces- of- section 

The motion of the test particle takes place on a three-dimensional subspace (defined by a 
particular value of C) embedded in the four-dimensional phase space, (X, X ,Y,Y). The usual 
two-dimensional surface-of-section (s-o-s) is defined by 

Y = 0, Y > 0, (5) 

and the coordinates on the section are (X,X). The geometrical interpretation is straightforward: 
we plot the X coordinate and velocity of the test particle at every conjunction with Neptune. For 
orbits exterior to Neptune (strictly speaking, with mean motion smaller than Neptune's), this 
occurs every time the test particle is aligned with the Sun and Neptune and is on the opposite 
side of the Sun from Neptune; for orbits interior to Neptune (i.e. with mean motion larger than 
Neptune's), this occurs every time the test particle is aligned with the Sun and Neptune and is on 
the same side of the Sun as Neptune. For the Solar system beyond Neptune, only the exterior 
orbits are of interest. 

In the s-o-s so defined, periodic orbits of the test particle appear as discrete points. The 
successive crossings of this surface by a quasiperiodic orbit lie on a closed smooth curve, while a 
chaotic orbit fills a two-dimensional area. 

As the regions in the neighborhood of orbital resonances are of particular interest here, the 
following remarks are in order. 

• The location of a (j + k) : j orbital period resonance is defined by Kepler's relation between 
the orbital period and the semimajor axis: 

2/3 

(6) 

where j and k are integers. 

• In general, at any resonance, one finds a family of stable periodic orbits that can be 
parametrized by the orbital eccentricity (or, equivalently, the Jacobi integral, if we use 
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Eqn. |6] in Eqn. |J) such that the semimajor axis equals a res for all orbits in this family, and 
the eccentricity ranges from zero to some maximum value, e max ; for e > e max the periodic 
orbit becomes unstable. 

• In some neighborhood, ±Aa of dyes 5 iiGcir each resonance, there exist stable quasiperiodic 
orbits that librate with finite amplitude about the exact resonant orbit; the half- width Aa 
depends upon the mass parameter /i, the resonance itself (i.e. a res ) and the Jacobi integral 
C. 



Let us consider the Neptune-Pluto 3:2 resonance as an illustrative example. If Pluto's 
inclination is neglected, and if one assumes a fixed circular orbit for Neptune, then the Jacobi 
integral C = 2.9798 for the Sun+Neptune+Pluto restricted three-body problem. Figure 1(a) 
is a s-o-s in the neighborhood of the 3:2 Neptune resonance for this value of C. This section 
shows examples of all three types of orbits. The center of the smooth curves lying near 
(X, X) = (—1.65, 0) corresponds to the periodic orbit whose period is exactly 3/2 that of Neptune. 
The smooth curves are quasiperiodic orbits surrounding this exact periodic orbit and represent 
test particle orbits that librate about the exact resonance; the size of a smooth curve bears a 
direct relationship to the amplitude of libration. Such librating orbits are phase-protected from 
having close encounters with Neptune even if the orbit is Neptune-crossing. It can be seen that, 
in general, the regions of quasiperiodic orbits are surrounded by chaotic orbits, i.e. beyond a 
certain libration amplitude, the smooth curves dissolve into a chaotic zone. Embedded within the 
"chaotic sea" are other stable zones or chains-of-islands that represent other nearby higher order 
and secondary resonances. The origin of the chaos is due to the overlap of the higher-order and 
secondary resonances. 

It is perhaps more meaningful to visualize the resonance zone in variables that can be more 
directly translated into "resonance width". The latter is traditionally (and somewhat loosely) 
defined as the range in semimajor axis where the orbital perturbations are large. In the restricted 
three-body problem with ^<1, this is the range in which the amplitude of the orbital parameter 
variations of the test particle is not simply linearly proportional to fj,, but is 0([i u ) with v ~ 0.5. 
In general, in the neighborhood of any resonance, we find regions of stable, librating orbits 
surrounded by a substantial chaotic zone where test particles may suffer close encounters with 
Neptune. In this paper, I will define the resonance half-width Aa as the maximum amplitude of 
the semimajor axis variations of stably librating resonant orbits. The value of Aa will be estimated 
from the s-o-s as described below. 

The phase space structure seen in the (X, X) plane can also be seen in other (generalized) 
canonical coordinates. For example, in plane polar coordinates the equivalent s-o-s (r,r,9,9) 
in the inertial reference frame is (r, r) with the same section condition as in Eqn. [|. The s-o-s 
condition can be rewritten in terms of polar coordinates as follows: 



0'=l " , B<1 




-7- 



where 9' is the longitude of Neptune, and n is the mean motion of the test particle. In other 
words, the (r, r) plane is equivalent to the (X,X) plane with points plotted at every conjunction 
of the test particle with Neptune. The (r, r) s-o-s is shown in Figure 1(b). 

Carrying this one step further by means of a canonical transformation to Delaunay 
variables, (M, J, A, L) where M, A are the mean anomaly and mean longitude and 
J = \/a(l — v 7 ! — e 2 ),L = \fa{l — e' 2 ) (cf. chapter XVII in Brouwer & Clemence (1961) 



and chapter 10 in Goldstein (1980) ), this same section can also be seen in the (M,J) plane. 
The (M, J) s-o-s is shown in Figure 1(c). It is evident from the latter figure that for Pluto-like 
resonance-locked orbits (i.e. for the particular value of C chosen in this s-o-s) the mean anomaly 
at conjunctions with Neptune librates about 180° and its maximum libration amplitude is ~100°. 
Orbits with larger amplitudes are chaotic. 

In addition to the maximum amplitude of M, the resonance libration region is bounded in 
J as well. Now, if we note that (i) the Jacobi integral provides a relationship between a and e 
(see Eqn. (|j)), and (ii) for a particular orbital resonance defined by the ratio (j + k) : j of the 
orbital periods of the test particle and Neptune, the librating orbits all have a phase-averaged 
semimajor axis (a) = a res , then any s-o-s can be labeled uniquely either by the value of C or, in 
the neighborhood of a particular orbital resonance, by the value of (e) (with the understanding 
that (a) = a res ). For example, for the s-o-s in Figure 1, we have C = 2.9798; near the 3:2 
resonance, (a) = 1.3104; therefore, from Eqn. @) we have (e) = 0.250. Thus, the phase space 
structure in the neighborhood of a particular orbital resonance can be systematically studied in 
a set of surfaces-of-section where all of them have a common value of (a) but various values of 
(e) (hence different values of C). This allows a visualization of the phase space near an orbital 
resonance as a function of the mean orbital eccentricity. 

Finally, in Figure 1(d) I have plotted the values of a against M for the same points 
represented in Figures l(a-c). Note that this is not a surface- of- section. One might call it a 
"pseudo-surface-of-section" , for a is not a canonical variable. However, the constancy of C imposes 
a relationship between a and e (cf. Eqn. ^|) which ensures that this plot looks very similar to the 
(M, J) s-o-s in Figure 1 (b) and the stable libration regions and the chaotic regions are easily 
distinguished. The reason for choosing this non-canonical variable is that the resonance width 
Aa of the stable libration zone can be readily determined from this figure. Therefore, the (a, M) 
plane — the pseudo-surface-of-section — will be used in the next section for a systematic look at 
the properties of the phase space in the neighborhood of several Neptune resonances. 

The libration timescale in the libration zone and the timescale for orbital instability in the 
chaotic zone are of particular interest for the dynamics of Kuiper Belt objects. These libration 
timescales have been determined for individual resonances and are also given in the next section. 
Typically, the timescale for orbital instability inside the chaotic zone surrounding a resonance is 
only a few times the libration period. However, close to the boundary between the libration zone 
and the chaotic zone, this timescale can become exceedingly large. In particular, for Pluto-like 
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orbits with eccentricity ~ 0.25 in the 3:2 Neptune resonance, the timescale for instability in the 
chaotic zone (i.e. for amplitude of the resonance angle greater than ~ 140°) is 0(1O 5 ) yr, but close 
to the edge of the libration zone this timescale increases by an order of magnitude or more. This 
is illustrated in Figure 2. 



3. Neptune Resonance Zones 



In principle, there is an infinite number of orbital resonances in the three-body problem. 
However, we will see that only a handful warrant any detailed investigation. 

Low-eccentricity orbits in the immediate neighborhood of a planet are unstable. There is 
a simple explanation for this orbital instability that follows from Chirikov's (1979) "resonance 
overlap criterion": first order (j + 1) : j resonances with (j + 1) > 0.51/z~ 2 / 7 overlap completely 
IjWisdom (1980) | . [Numerical experiments by Duncan et al. (1989) show that the numerical 



coefficient on the right hand side in this relation must be revised to 0.44]. Complete overlap of 
neighboring resonances results in the destabilization of all periodic and quasi-periodic orbits in 
the neighborhood of those resonances. In other words, circular test particle orbits are unstable in 
the immediate neighborhood \a — 1| ^ 1.5/U 2 / 7 of a planet's orbit. For Neptune (/i = 5.146 x 10~ 5 ), 
this criterion shows that all first order resonances with (J + 1) > 8 are completely overlapping, so 
that circular orbits beyond Neptune with orbital radii less than about 33 AU are unstable. 

Therefore, beyond the semimajor axis range of about 33 AU, at most only seven first order 
(j + 1) : j resonances with j = 1, 2, ...7 are isolated from each other (at least for low eccentricities). 
However, I have found numerically that low-eccentricity orbits at the 8:7 and 7:6 Neptune 
resonances are also chaotic. This is not entirely surprising, for the resonance overlap criterion is 
quite approximate, and Wisdom's scaling law is strictly valid only in the asymptotic limit j>l. 
Therefore, only five first order (j + 1) : j resonances with j = 1, 2, ...5 are of potential significance 
for the long-term storage of Kuiper Belt objects. 

Furthermore, since all second and higher order resonances in the first-order resonance overlap 
region are also destroyed, it follows that there are only seven second order Neptune resonances of 
potential interest: (J + 2) : j with j = 1, 3, 5, 13. Of these, I have found numerically that only 
three (3:1, 5:3, and 7:5) have resonance widths Aa J> 0.005. Therefore, the phase space in the 
neighborhood of only these three second order resonances will be discussed here. 

Figures 3-10 show the properties of these first and second order Neptune resonances. They 
are arranged in order of decreasing mean heliocentric distance. Each figure has seven parts: the 
first four [(a)-(d)] show the resonance zones in the pseudo-surfaces-of-section [plots of (a, M)] for 
(e) = 0.1,0.2,0.3 and 0.4, respectively; note that the vertical scale is the same in all these figures, 
so that the variations in the extent of the libration and chaotic zones with (e) and with the mean 
distance of the resonance from Neptune can be readily discerned. The next two panels [(e) and 
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(f)] show examples of librating, resonance-locked orbits in the rotating reference frame. The last 
one, panel (g), shows the period of libration as a function of the mean orbital eccentricity and the 
amplitude of libration of the resonance angle, (ft. A few words are in order regarding the latter 
variable. 

For a (j + k) : j resonance the resonance angle is defined as follows: 

cj)=(j + k)X - jX' - km (7) 

where A and w are the mean longitude and longitude of perihelion of the test particle, and A' 
is the mean longitude of Neptune. (Note that j is a positive integer; A: is a negative integer 
for interior resonances, positive for exterior resonances; the latter are the only resonances of 
interest in this paper.) (ft is the natural variable that arises in the perturbative analysis of orbital 
resonances, and it can be used to make a first approximation theory that models the resonant 
motion of the test particle with a single degree-of-freedom pendulum-like dynamical system [see 
e.g. Malhotra (1994)| . For a test particle locked in a stable exterior resonance, (ft librates about 



a mean value which is usually 180°. [There are two exceptions to this: (i) for sufficiently high 
eccentricity, librations of the resonance angle about 0° are also possible; and (ii) the 2:1 and the 
3:1 resonances allow for asymmetric librations where the center of libration of the resonance angle 
is displaced away from 180°. These are discussed further below.] At conjunctions with Neptune, 
the resonance angle is related to the mean anomaly: <f> = jM + 0(e). Therefore, the librations or 
chaotic variations of M evident in the surfaces-of-section are also reflected in the behavior of (ft. 

The physical significance of the resonance angle can be seen by noting that when the test 
particle is at perihelion, A = zu, so that 

(ft = j(w — A') (at perihelion). 

Thus, the behavior of (ft gives a direct measure of the longitude separation of Neptune from the 
test particle's perihelion. This is obviously a critical quantity, especially for high-eccentricity 
Neptune-crossing orbits, as the stable libration of (ft then ensures that the particle is protected 
from catastrophic close encounters with Neptune. The maximum excursion from the center of 
libration of (ft then gives a measure of the smallest possible separation between Neptune and the 
test particle's perihelion in stable resonance-protected orbits. 

The resonance libration period as a function of the amplitude of libration of (ft is shown for 
various values of (e) in the last panel (g) of Figures 3-10. These curves are rather counter-intuitive. 
Recall that in the oft-used pendulum model for non-linear resonances, the pendulum libration 
period increases monotonically with amplitude and becomes arbitrarily large close to the separatrix 
(i.e. the orbit that separates the oscillations from rotations of the pendulum). In contrast, in 
Figures 3-10(g), we see that, except for the asymmetric librations at the 2:1 and 3:1 resonances, 
the libration period in general decreases with increasing libration amplitude. Furthermore, near a 
k-th order resonance, the leading order (in eccentricity) resonance term in the perturbing potential 
for the test particle is proportional to the k-th power of the eccentricity, hence the small-amplitude 
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libration period might be expected to vary as ~ e~ fc / 2 ; this analogy from the pendulum model 
does not hold either, as evident in Figures 3-10(g). Thus, the standard non-linear pendulum 
model should be used with caution for orbital resonances. 

The general characteristics of the phase space near these first and second order Neptune 
resonances are summarized as follows. (Note that the 2:1 and the 3:1 resonances are exceptional 
in many respects, and are discussed later in detail.) 

[cf. Figures 3-10 (a-d)] In general, the phase space in the neighborhood of a resonance 
contains a zone of stably librating, resonance-locked orbits surrounded by a zone of chaotic 
orbits. The width of the resonance zone (including the stable as well as the chaotic regions) 
increases with (e); however, the width of the libration zone, in general, shrinks as the 
resonance separatrix dissolves into a chaotic layer that increases in thickness with (e). The 
typical width of a Neptune resonance libration zone is 2Aa « 0.02, or about 0.6 AU. 

[cf. Figures 3-10 (e)] The periodic orbit at exact (j + k) : j resonance has a j-fold symmetry 
in the rotating frame: it makes j perihelion passages during a synodic period. The perihelion 
(and aphelion) longitudes are spaced 360°/j. The quasi-periodic orbits librating about the 
exact resonance also exhibit j-fold symmetry in the rotating frame when traced over a 
complete libration period. For the stable librators (i.e. resonance angle <j) librating about 
180°), the opposition with Neptune occurs near the test particle's aphelion; hence these may 
be called "aphelion librators". 

[cf. Figures 6,7,9,10 (f)] For sufficiently large (e), a new type of periodic (resonant) orbit 
is possible, where one of the j perihelion passages during a synodic period occurs near 
Neptune. In fact, the perihelion loop encompasses Neptune. For this orbit and for the 
quasi-periodic orbits in its neighborhood, the resonance angle <fi librates about 0°. These 
orbits may be called "perihelion librators" . However, these perihelion librators are probably 
not of practical importance, as unmodeled perturbations are likely to destabilize them. 

[cf. Figures 3-10 (g)] In general, the libration periods are smaller for resonances located 
closer to Neptune, and the libration periods decrease with libration amplitude. Typical 
small amplitude libration periods are (1-2) xlO 4 yr. 

The exceptional cases of the 2:1 and the 3:1 resonances are discussed below. 

The 2:1 Neptune resonance 

The synodic period of a test particle near the 2:1 resonance is approximately equal to its 
orbital period. Thus, conjunctions with Neptune occur every ~ 330 years. The properties of the 
phase space near the 2:1 Neptune resonance are summarized as follows (cf. Figure 4). 
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The resonance half-width, Aa, increases — from ~ 0.001 to ~ 0.018 — as (e) increases from 
~ to 0.4. The phase space at this resonance is largely regular; there is no discernible 
chaotic zone until the eccentricity exceeds ~ 0.25. 

For (e) exceeding a critical value (which is ~ 0.04 for the 2:1 Neptune resonance), there 
are two types of librating orbits: (i) symmetric librators whose perihelia librate (with large 
amplitude) with a mean value 180° away from Neptune, and (ii) asymmetric librators 
whose perihelia librate (with amplitude smaller than ~ 90°) about a mean value which 
depends upon (e). (Surprisingly, the orbit with aphelion exactly at Neptune's longitude is 
unstable: it coincides with the separatrix that divides the phase space into the symmetric 
and asymmetric libration zones.) The existence of such asymmetric librations has been 



noted previously (e.g. |Message (1958)1 , [Beauge (1994)| ), but the work here provides the first 



quantitative analysis for the specific case of Neptune's resonances. 

Figures 4(e,f) show examples of these two types of orbits. Note that for the asymmetric 
librators, there are two independent centers of libration, offset equally on either side of 
Neptune's longitude. Upon inspection of the behavior of the resonance angle </>, it is found 
that these centers of perihelion libration (relative to Neptune) vary from (±)180° for (e) ~ 
to (±)50° for (e) « 0.4. 

The phase space for the symmetric librators shrinks rapidly with increasing (e), while that 
for the asymmetric librators increases slightly. 

At small values of (e), there is no evidence for a chaotic zone near this resonance, but for 
(e) > 0.25, the largest amplitude librators become chaotic. A thin chaotic zone also appears 
in the neighborhood of the separatrix that separates the asymmetric librations from the 
symmetric librations. 

The eccentricity variations of the resonance-locked orbits have amplitude be < 0.01. (Indeed, 
within the planar circular restricted three-body model, an initially circular test particle 
orbit at the 2:1 resonance has its eccentricity pumped up to only ~ 0.03.) The eccentricity 
variations on the chaotic orbits are generally of considerably larger magnitude. 

Figure 4(g) shows libration amplitude vs. libration period. For the asymmetric (small 
amplitude) librators, the libration period decreases monotonically with (e) and increases 
with libration amplitude (reminiscent of the standard non-linear pendulum). The libration 
period becomes indefinitely large at the separatrix between the asymmetric and symmetric 
librators. The small amplitude libration periods are several tens-of-thousands of years long. 



The 3:1 Neptune resonance 



The synodic period of a test particle near the 3:1 resonance is half its orbital period. Thus, 
conjunctions with Neptune occur approximately every 248 years. The properties of the phase 
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space near the 3:1 Neptune resonance (cf. Figure 3) are very similar to that near the 2:1 resonance; 
some specifics are noted below. 

- The resonance half-width, Aa, increases — from < 0.001 to ~ 0.016 — as (e) increases from 
~ to 0.4. The phase space at this resonance is regular; there is no discernible chaotic zone 
for this range of eccentricities. 

- For librating orbits, perihelion and aphelion occur at approximately the same longitude 
relative to Neptune. As in the 2:1 resonance, there are two types of librating orbits - the 
symmetric and asymmetric librators; however, the asymmetric librations appear only at (e) 
exceeding ~ 0.13. 

- The small amplitude libration periods are very long (~ 0.4 Myr) for low eccentricity e ^ 0.1 
orbits, but considerably shorter (30,000-80,000 yr) for higher eccentricities. 



4. Summary and Discussion 



Numerical studies of the stability of low eccentricity, low inclination orbits of small objects in 
the trans-Neptunian Kuiper Belt subject to the gravitational perturbations of the giant planets 
have shown that the inner edge of the Kuiper Belt is at about 34 AU heliocentric distance; beyond 
34 AU and up to about 42 AU there are interspersed regions of stability and instability [Levison & 



Duncan (1993) , Holman fc Wisdom (1993) 1- This structure bears a complex correlation with the 



locations of Neptune orbital resonances. Other studies | Malhotra (1993)1 , Malhotra (1995) | suggest 
that early in the history of the Solar system, the majority of Kuiper Belt objects were swept into 
eccentric orbits in narrow zones located at Neptune's orbital resonances, and that the regions 
in-between the resonances would have been largely cleared of residual planetesimals. Recent and 
ongoing observational surveys of the outer Solar system indicate the presence of a large population 
of small bodies beyond Neptune [ Jewitt fc Luu (1995) | . These are quite likely the source of short 
period comets [Duncan et al. (1988), Quinn et al. (1990)]. Furthermore, their orbital distribution 
is likely to hold clues to the formation and early dynamical evolution of the outer Solar system 
jMalhotra (1995)| |. All of these considerat ions have motivated the present study. 

Using the planar circular restricted three body model (with the Sun, Neptune and a test 
particle), I have described in this paper the basic phase space structure in the neighborhood of 
Neptune's exterior orbital resonances. The details may be gleaned from Figures 3-10. A succinct 
summary of these results is given in Figure 11 which shows the locations and widths of the stable 
resonance libration zones. In general, these stable zones are bounded by chaotic layers of thickness 
that generally increases with eccentricity and decreases with distance from Neptune. (As the mean 
semimajor axis and eccentricity of chaotic orbits is not well defined, it is not possible to represent 
the chaotic layers in such a figure.) 
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The planar circular restricted three body model is the simplest dynamical model for the orbital 
dynamics of small objects in the Kuiper Belt. Although this model may appear oversimplified, it 
provides a reasonable description and explanation for much of the dynamical behaviors found in 
the numerical studies mentioned above.[] For example, the location of the inner edge of this Belt 
at ~ (33-34) AU is readily understood in terms of orbital instability induced by overlapping first 
order Neptune resonances; the additional perturbations due to Neptune's non-circular, inclined, 
and time-varying orbit and the perturbations of the other planets do not significantly change the 
location of this boundary. 

The relatively isolated first and second order Neptune resonances beyond 34 AU (in semimajor 
axis) provide narrow stable libration regions for the long term storage of Kuiper Belt objects in 
eccentric (often Neptune-crossing) orbits. This paper provides first approximation estimates for 
the locations and widths of these regions, and for the dynamical properties of resonant orbits. 
It also provides a direct visualization — in two-dimensional surfaces-of-section - of the global 
phase space structure (i.e. the chaotic and stable regions) in the vicinity of orbital resonances A 
characteristic of stable resonant orbits is the libration of the perihelion about a longitude well 
removed from Neptune's location. Typical libration periods are several tens-of-thousands of years. 
The libration zones are generally surrounded by narrow chaotic zones where orbits are unstable on 
timescales of a few libration periods, or ~ 10 5 yr. At the boundary between the stable resonance 
libration zone and the chaotic zone, the timescale for orbital instability may be exceedingly long. 
This paper also provides a quantitative analysis of libration periods and their dependence on 
libration amplitudes and orbital eccentricities at all the major exterior mean motion resonances of 
Neptune. 

Of course, the planar circular restricted three body model does not provide a complete 
picture. Indeed, I expect that it underestimates the extent of the chaotic zones and overestimates 
the sizes of the stable libration zones near Neptune resonances. But perhaps the most important 
missing element is the effect of secular resonances [see e.g. Knezevic et al. (1991)] on the long 



term dynamics of Kuiper Belt objects. For example, in some regions, secular resonances produce 
significant inclination excitation; in such cases, the assumption of planarity becomes a poor 
approximation. [Pluto's orbit is a case in point ( Nacozy fc Diehl (1978j ).] In order to discern these 



effects, one has to build much more elaborate analytical models or resort to extensive numerical 
integrations. 



As this work was being completed, I became aware of two preprints | Morbidelli et al. (1995) 



and Duncan et al. (1995) | on the same subject. These use different approaches to the problem: 



the present work has focussed on the phase space near mean motion resonances with Neptune only 



and gives results that are exact but for a highly simplified physical model; Morbidelli et al. (1995) 



1 Those numerical studies assumed a realistic three-dimensional physical model in which test particles were 
perturbed by all the four giant planets, Jupiter, Saturn, Uranus and Neptune; the planets' orbits were integrated 
self-consistently under their mutual gravitational perturbations. 
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discuss the dynamics near mean motion resonances as well as secular resonances using primarily 
approximate semi-analytic models; and Duncan et al. (1995) discuss the dynamics near mean 
motion resonances, secular resonances as well as non-resonant regions using a purely numerical 
approach. A comparison of the present work with these two preprints follows. 

Morbidelli et al. (1995) describe semi-analytic and numerical investigations that include the 
effects of all four giant planets on the dynamics of test particles in the Kuiper Belt. A comparison 
of their results for the locations and widths of mean motion resonances (see their Figure 1) with 
the estimates obtained here (cf. Figure 11) shows that their models significantly overestimate 
the sizes of the resonance libration regions. Morbidelli et al.'s estimates of the resonance widths 
were obtained by averaging the planar, circular restricted problem of the four giant planets plus 
a test particle over fast variables, retaining only the dependence on the resonance angle. Thus, 
their analysis does not account for the overlapping secondary resonances at the edges of the mean 
motion resonance libration zones that produce a chaotic layer and cause the libration regions to 
shrink rapidly with increasing eccentricity. This is the explanation for the differences between 
their figure 1 and figure 11 here. In addition, Morbidelli et al. find a singularity in the averaged 
Hamiltonian that causes a great increase in the resonance width near Neptune-crossing values of 
the eccentricity. However, this is an artifact, and does not appear in the unaveraged problem, 
as is obvious from the surfaces-of-section shown in the present work. The technique I used here 
does not make any averaging approximations, but uses the full unaveraged Hamiltonian for the 
circular planar restricted problem, albeit without the secular effects of Jupiter, Saturn and Uranus. 
Morbidelli et aVs models also show that secular effects are negligible near mean motion resonances 
except in the case of the 3:2 Neptune resonance. In the special case of the 3:2 resonance, as 
mentioned above, Knezevic et al. (1991) had previously shown the existence of the z^is secular 
resonance embedded within the mean motion resonance; in addition, analysis of Pluto's orbit 
had also previously revealed the existence of another resonance characterized by the libration of 
the argument of perihelion flNacozy &: Diehl (1978)| ). Both of these resonances have the effect 
of exciting the inclination of the test particle. These resonances are not modeled in the present 
work, but Morbidelli et al. describe in some detail the location and widths of both these and other 
secular resonances. However, note that because most regions near secular resonances are regions 
of orbital instability, they are not likely to be of importance as reservoirs of Kuiper Belt objects 
although they may drive their transport to small heliocentric distances at the present epoch. 

Duncan et al. (1995) have performed quite extensive numerical integrations of test particle 
orbits subject to the gravitational perturbations of the four giant planets. They describe their 
results in terms of the stability (over their integration time of 1-4 Byr) of orbits with particular 
initial semimajor axes, eccentricities and inclinations. Their main conclusions are consistent with 
the results of the present work: namely, that particles with low initial inclination (~ 1°) and 
initial perihelion distance less than ~ 35 AU are unstable on timescales short compared to the 
age of the Solar system, except that particles librating in low order mean motion resonances with 
Neptune remain phase-protected from close encounters with that planet. They have also analyzed 
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the 3:2 Neptune resonance in more detail for the particular case with initial eccentricities of 0.2 
and found three regimes: stable orbits deep in resonance with small libration amplitudes, < 70°; 
an intermediate region with libration amplitudes in the range 70° — 130° where the timescale 
for instability is ~ 10 9 yr; and highly unstable orbits for libration amplitudes exceeding ~ 130° 
with stability timescales of < 10 8 yr. Similar dynamics is found for all the Neptune resonances 
studied in this paper. In comparison, the present work has shown that the region of stability 
deep in resonance as well as the highly chaotic region at large amplitudes are both explained by 
the perturbations of Neptune alone (assumed to be on a circular orbit); the intermediate regime 
where the timescale for instability is ~ 10 9 yr is very small in the restricted three body model but 
evidently expands greatly (at the expense of the stable libration zone) in the full N-body model 
for the outer Solar system. [In this context, we note that if Malhotra's (1995) theory for the 
orbital distribution of Kuiper Belt objects is correct, then this intermediate libration amplitude 
regime which allows a long-term leakage of resonant objects may be the primary source of the 
Jupiter-family short-period comets.] Duncan et al. also find that the resonance-protection fails 
only at inclinations exceeding ~ 25°. Thus, we conclude that the planar circular restricted three 
body model used in the present work provides a fairly good description of these main trends. A 
clear advantage of this simple model is that its analysis is quite "inexpensive" compared to the 
months of CPU time expended in the numerical simulations by Duncan et al. Furthermore, the 
numerical integration of orbits for the age of the Solar system at sufficiently high resolution in the 
space of initial conditions remains beyond the reach of present-day computers, as acknowledged 
by Duncan et al. Thus, the global picture of the phase space structure obtained in the present 
work provides an inexpensive yet reasonably good baseline for understanding the dynamics in the 
trans-Neptunian Solar system, and is particularly useful near mean motion resonances. 

This research was done while the author was a Staff Scientist at the Lunar and Planetary 
Institute which is operated by the Universities Space Research Association under contract 
no. NASW-4574 with the National Aeronautics and Space Administration. Partial support for 
this work was also provided by NASA's Origins of Solar Systems Research Program under grant 
no. NAGW-4474. This paper is Lunar and Planetary Institute Contribution no. 870. 
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FIGURE CAPTIONS 

Fig. 1 — The phase space structure in the neighborhood of the Neptune-Pluto 3:2 orbital resonance 
determined from the circular planar restricted three-body model. The same surface-of-section 
is shown in different variables: (a) cartesian variables, (X, X); (b) plane polar variables, (r, f); 
(c) Delaunay variables (mean anomaly M, canonical momentum J = \/a(\ — \/l — e 2 )). (d) is 
a pseudo-surface-of-section since the semimajor axis, a, is not a canonical variable. However, 
qualitatively it is very similar to the surface-of-section (c), and it provides a visualization of the 
resonance width in terms of the semimajor axis and the amplitude of the perihelion libration. 

Fig. 2 — The time variation of the semimajor axis for orbits with e ~ 0.25 near the boundary 
between the stable libration zone and the chaotic zone at the 3:2 Neptune resonance (cf. Figure 1). 
(The unit of length is Neptune's semimajor axis, and the unit of time is Neptune's orbital period, 
approximately 165 yr.) The Jacobi integral for all four orbits is C = 2.9798; the initial conditions 
(X,X,Y,Y) for the four orbits are (a)— (-1.6100,0.0,0.0,0.92442), (b)— (-1.6105,0.0,0.0,0.92508), 

(c) — (-1.6106,0.0,0.0,0.92521), (d)— (-1.6110,0.0,0.0,0.92574). These orbits are all chaotic and are 
shown in order of decreasing distance to the stability boundary. The orbits in (a), (b) and (c) 
become unstable after 7, 31 and ~ 200 librations (i.e. ~ 8 x 10 4 , ~ 3 x 10 5 and ~ 2 x 10 6 yr), 
respectively. 

Fig. 3 — The 3:1 Neptune resonance located at a res = 2.08008 (i.e. 62.6 AU). The phase space 
for orbits with (a) C = 3.3508, (e) = 0.1; (b) C = 3.3070, (e) = 0.2; (c) C = 3.2324, (e) = 0.3; 

(d) C = 3.1244, (e) = 0.4. The phase space is largely regular at this resonance. The width of 
the libration zone increases in with (e). For (e) ^> 0.13, the libration zone splits into two zones 
of asymmetric librations surrounded by a narrow zone of large-amplitude symmetric librations. 
Panel (e) shows a symmetrically librating orbit (with e « 0.1) in the rotating frame; panel (f) 
shows an asymmetrically librating orbit (with e « 0.2). The symbols and + indicate the fixed 
locations of the Sun and Neptune in the rotating frame. Panel (g) shows the libration period (in 
units of Neptune's orbital period) as a function of the maximum excursion of the resonance angle 
(f) from 180°. Note that the center of the asymmetric librators moves away from "</> = 180°" with 
increasing (e), as indicated by the beginning of the curves labeled by the different values of (e). 
The libration period increases with libration amplitude for the asymmetric librators, but decreases 
with amplitude for the symmetric (large amplitude) librators; at the separatrix on the boundary 
between the symmetric and asymmetric librations, the libration period becomes indefinitely large 
(although this is not quite obvious due to the finite resolution in this figure). 

Fig. 4 — The 2:1 Neptune resonance located at a res = 1.5873 (i.e. 47.8 AU). The phase space 
for orbits with (a) C = 3.1272, (e) = 0.1; (b) C = 3.0989, (e) = 0.2; (c) C = 3.0338, (e) = 0.3; 
(d) C = 2.9395, (e) = 0.4. The phase space is largely regular at this resonance; a chaotic zone 
appears at large amplitude librations for eccentricity exceeding ~ 0.25. The width of the libration 
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zone increases in with (e). For (e) > 0.04, the libration zone splits into two zones of asymmetric 
librations surrounded by a narrow zone of large-amplitude symmetric librations. Panels (e) and 
(f) show a small-amplitude asymmetrically librating orbit and a large amplitude symmetrically 
librating orbit (both with e ~ 0.3) in the rotating frame. (The motion of the test particle over 
a complete libration period is traced in these figures.) Note that in both types of orbits the 
perihelion occurs away from Neptune's longitude. Panel (g) shows the libration period (in units of 
Neptune's orbital period) as a function of the maximum excursion of the resonance angle <p from 
180°. The center of the asymmetric librators is indicated by the beginning of the curves labeled 
by the different values of (e). The libration period increases with libration amplitude for the 
asymmetric librators, but decreases with amplitude for the symmetric (large amplitude) librators; 
at the separatrix on the boundary between the symmetric and asymmetric librations, the libration 
period becomes indefinitely large. 

Fig. 5 — The 5:3 Neptune resonance located at a res = 1.4057 (i.e. 42.3 AU). The phase space 
for orbits with (a) C = 3.0709, (e) = 0.1; (b) C = 3.0349, (e) = 0.2; (c) C = 2.9739, (e) = 0.3; 
(d) C = 2.8848, (e) = 0.4. The width of the libration zone increases for (e) up to ~ 0.2, but 
then shrinks with increasing (e) as the larger amplitude resonant orbits dissolve into the chaotic 
zone. Note that for (e) > 0.3, the formerly period-2 exact resonant orbit in the surface-of-section 
becomes a period-3 orbit; this is merely a reflection of the fact that at large, Neptune-crossing 
eccentricities the perihelion loop results in two conjunctions with Neptune in rapid succession 
[cf. panel (c) and (f)]. Panels (c) and (f) show the periodic orbit at exact resonance, and an 
orbit librating with finite amplitude, respectively (both with e ~ 0.3) in the rotating frame. (The 
motion of the test particle over a complete libration period is traced in these figures.) Note that 
in both types of orbits the perihelion occurs away from Neptune's longitude. Panel (g) shows the 
libration period (in units of Neptune's orbital period) as a function of the libration amplitude of 
the resonance angle <f>. The libration period increases with libration amplitude only for small (e); 
it decreases with amplitude for larger (e). 

Fig. 6 — The 3:2 Neptune resonance located at a res = 1.3104 (i.e. 39.4 AU). The phase space 
for orbits with (a) C = 3.0411, (e) = 0.1; (b) C = 3.0065, (e) = 0.2; (c) C = 2.9470, (e) = 0.3; 
(d) C = 2.8616, (e) = 0.4. The width of the main libration zone decreases with increasing (e), 
as the chaotic zone surrounding the libration region expands. Panel (e) shows an example of 
a stably librating orbit (with e ~ 0.3) in the rotating frame; note that this is an "aphelion 
librator": the test particle is near aphelion when passing Neptune's longitude, and its perihelion 
librates about ±90° away from Neptune. For (e) > 0.3, a new libration zone appears in which 
the librating orbits have a perihelion near Neptune. This is indicated by the appearance of a 
chain of two libration islands (distinct from the main resonance zone) in panels (c)-(d); these are 
the "perihelion librators" [cf. panel(f)]; The new libration zone increases with increasing (e). An 
example of a perihelion librator is shown in panel (f); note that one of the two perihelion loops 
in each synodic period encompasses Neptune. Panel (g) shows the libration period (in units of 
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Neptune's orbital period) of the aphelion librators as a function of the libration amplitude of the 
resonance angle <p. The libration period generally decreases with libration amplitude, except for 
very small eccentricity (e < 0.1), large-amplitude librators. 

Fig. 7 — The 7:5 Neptune resonance located at a res = 1.2515 {i.e. 37.7 AU). The phase space for 
orbits with (a) C = 3.0255, (e) = 0.1; (b) C = 2.9913, (e) = 0.2; (c) C = 2.9334, (e) = 0.3; (d) 
C = 2.8496, (e) = 0.4. The width of the main libration zone (for aphelion librators) decreases 
rapidly with increasing (e), as the chaotic zone surrounding the libration region expands. In fact, 
for (e) > 0.3 [cf. panel (d)] the aphelion libration zone has completely disappeared. A libration 
zone for perihelion librators appears at e ~ 0.3 and expands slightly with increasing (e). Panels 
(e) and (f) show examples of an aphelion librator with e ~ 0.2, and a perihelion librator with 
e ps 0.3, respectively. Panel (g) shows the libration period (in units of Neptune's orbital period) of 
aphelion librators as a function of the libration amplitude of the resonance angle 4>. The libration 
period generally decreases with libration amplitude, except for very small eccentricity (e < 0.1), 
large-ampitude librators. 

Fig. 8 — The neighborhood of the Neptune 4:3 orbital resonance located at a res = 1.2114 (i.e. 36.5 
AU). The phase space for orbits with (a) C = 3.0158, (e) =0.1; (b) C = 2.9824, (e) = 0.2; (c) 
C = 2.9253, (e) = 0.3; (d) C = 2.8430, (e) = 0.4. The libration zone shrinks rapidly with (e) as the 
larger amplitude librators dissolve into the chaotic zone. For (e) > 0.175, the formerly period-1 
orbit corresponding to exact 4:3 resonance bifurcates into a period-2 orbit in the surface-of-section; 
this is merely a reflection of the fact that for a Neptune-crossing value of the eccentricity, the 
perihelion loop causes two conjunctions with Neptune in each synodic period. Panel (e) and (f) 
show examples of small amplitude librating orbit with e « 0.1 and e 0.2, respectively. Panel (g) 
shows the period of libration as a function of the libration amplitude of the resonance angle, 4> (for 
(e) = 0.1,0.2,0.3 and 0.4). The period as well as the maximum amplitude for stable librations 
decreases with increasing (e) . 

Fig. 9 — The neighborhood of the Neptune 5:4 orbital resonance located at a res = 1.1604. 
The phase space for orbits with (a) C = 3.0054, (e) = 0.1; (b) C = 2.9727, (e) = 0.2; (c) 
C = 2.9170, (e) = 0.3; (d) C = 2.8364, (e) = 0.4. Notice that the main libration zone (for aphelion 
librators) shrinks rapidly with (e). A libration zone for perihelion librators appears for (e) > 0.2, 
as indicated by the appearance of a chain of two libration islands in (c)-(d); this zone increases 
with increasing (e). Panel (e) shows the stable periodic orbit with e ~ 0.3 at the 5:4 resonance (an 
aphelion librator) in the rotating frame coordinates; panel (f) shows an example of a perihelion 
librator with e ~ 0.3. Notice the four-fold symmetry in the rotating frame. Panel (g) shows 
the period of libration for the aphelion librators as a function of the libration amplitude of the 
resonance angle (for (e) = 0.1,0.2,0.3 and 0.4) 

Fig. 10 — The neighborhood of the Neptune 6:5 orbital resonance located at a res = 1.1292 
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(i.e. 34.0 AU). The phase space for orbits with (a) C = 3.0002, (e) = 0.1; (b) C = 2.9680, (e) = 0.2; 
(c) C = 2.9131, (e) = 0.3; (d) C = 2.8334, (e) = 0.4. Notice that the main libration zone (for 
aphelion librators) shrinks rapidly with (e); for (e) > 0.1 the center of this libration zone bifurcates 
into a period-2 orbit; this is simply a reflection of the fact that the perihelion loop causes two 
conjunctions with Neptune in rapid succession during each synodic period. Also for (e) > 0.1, 
a new libration zone appears for perihelion librators [indicated by the appearance of a chain of 
two libration islands with centers at M = 0, it in (b)-(d)] whose extent increases with increasing 
(e). The increased complexity of this libration zone in (d) is a reflection of the fact that the 
perihelion loops for large amplitude librations are so large that there are four conjunctions with 
Neptune during each synodic period. Panels (e) and (f) show examples an aphelion librator and 
of a perihelion librator with e ps 0.3 are the 6:5 resonance. Notice the five-fold symmetry in the 
rotating frame. Panel (g) shows the period of libration as a function of the libration amplitude of 
the resonance angle, 4> (for (e) = 0.1,0.2,0.3 and 0.4) for the aphelion librators. 

Fig. 11 — The locations and widths of first and second order Neptune resonances in the Kuiper 
Belt, as determined by the planar, circular restricted three-body model. The shaded region on the 
extreme left indicates the chaotic zone of "resonance overlap" in the vicinity of Neptune's orbit. In 
the region above dotted line, orbits are Neptune-crossing; these orbits are dynamically short-lived 
(due to close encounters with Neptune), unless protected by orbital resonances. 
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